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Abstract 

The superconducting state of a two-dimensional d-p model is studied from the spin fluctuation point of 
view by using a strong coupling theory. The fluctuation exchange (FLEX) approximatoin is employed to 
calculate the spin fluctuations and the superconducting gap functions self-consistently in the optimal- and 
over-doped regions of hole concentration. The gap function has a symmetry of d x 2_ y 2 type and develops 
below the transition temperature T c more rapidly than in the BCS model. Its saturation value at the 
maximum is about IO^b^c- When the spin fluctuation-induced superconductivity is well stabilized at low 
temperatures in the optimal regime, the imaginary part of the antiferromagnetic spin susceptibility shows a 
very sharp resonance peak reminiscent of the 41 meV peak observed in the neutron scattering experiment 
on YBa2Cu307. The one-particle spectral density around k — (n,0) shows sharp quasi-particle peaks 
followed by dip and hump structures bearing resemblance to the features observed in the angle-resolved 
photoemission experiment. With increasing doping concentration these features gradually disappear. 

1 Introduction 

Among a number of mechanisms proposed for the high temperature superconductivity since its discovery in 
1986 the spin fluctuation mechanism has continued to be one of the most promising mechanisms, [jj], 0, |3|, |i| 
||, ^, Q Following systematic studies based on the parameterized spin fluctuation theory, more quantitative 
studies were carried out by using the Hubbard model §, |, [HJ 0, (0§ [H) and the d-p model |]| ^ 
within the so-called fluctuation exchange (FLEX) approximation. 

In a previous article [JTsJ (hereafter referred to as I) we have developed a strong coupling theory of spin 
fluctuation-induced superconductivity on a microscopic basis by using a two-dimensional d-p model and applied 
it to interpret high T c cuprates. We adopted the fluctuation exchange (FLEX) approximation in calculating 
the transition temperature T c , the dynamical susceptibility, nuclear spin-lattice relaxation rate, and the one- 
particle spectral density above T c . The results compared rather well with the experimental results on cuprates 
in the optimal- and the over-doped regimes. 

In the present article we extend the calculation to the superconducting state below T c . By using the d-p 
model with the same parameter values as in I, we solve the non-linear Dyson-Gor'kov equations numerically and 
calculate various physical quantities; the energy gap function, the one-particle density of states, dynamical 
susceptibility and the one-particle spectral density. The results of calculation are discussed in comparison 
with the existing experimental results, particularly the resonance peak at 41 meV, observed in the neutron 
scattering experiment on YBa 2 Cu30y [|l9[ |2(], ^l], [2j| and the observed angle-resolved photoemission 
spectrum (ARPES), showing a resolution- limited quasi-particle peak followed by a dip and a hump structure 
around the (tt,0) point. ||, || ||, |f| 

In what follows we first summarize the model and the formalism we use here in § 2. In § 3 the results of 
calculation are presented and discussed in comparison with the existing experimental results. Finally, § 4 is 
devoted for summary of conclusions and general discussion. 



* Japan Science and Technology Corporation, Domestic Research Fellow. 



1 



2 Model and Approach 



The model we use here is the same as in I; the d-p model consisting of the d x 2_ y 2 orbitals of copper atoms and 
the p a (p x and p y ) orbitals of oxygen atoms in the CuC>2 plane of cuprates. The energy levels of these orbitals 
are e d and e p , respectively and the transfer matrix elements are considered only between the neighboring d 
and p orbitals (t dp ) and between the neighboring p x and p y orbitals (t pp ) where the parities of both the d x 2_ y 2 
orbital and the p a (p x and p y ) orbitals are taken into account. For brevity the Hubbard type electron-electron 
interaction U is considered only within a d orbital, and the other effects are assumed to be contained in e d 
and e p as the Hartree potentials. The Hamiltonian may be written as follows: 



H — Hq + H\ 



(1) 



H = $><4*<fca+ J2[{e p -(-l) m t p (k)}pl k(7 p mka + {it m (k)4 a p mka + h.c.}]], (2) 

k.cr 



m—1 



with 



and 



t m (k) = V2t dp [sm(k x a/2) - (-l) m sm(k y a/2)}, 
t p (k) — At pp sin(fc x a/2) sm(k y a/2), 



(3) 



= ^EEE^U^-,!^!^, (4) 

q k k' 

where dka and p m ko{m =1,2) are the Fourier transforms of the annihilation operators for the electrons in the 
d and p m — -±=\p x — (— l) m p y ] orbitals, respectively, and No is the number of atoms in the crystal. 
We denote the diagonalized expression for H as follows: 
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Hq = ^2Y1 E ^ka^ mka a mka , (5) 



■ a m—1 



«lfccr 




011k 


012k 


013k 




dka 


l2ka 




021k 


022k 


023k 




Piker 


»3fccr 




031k 


032k 


033k 




_ P2ka 



Then the Green's functions for the non-interacting system is given by 
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G { V{k,iu n ) = m \k0muk 1 , 



(6) 



(7) 



where A(V) = 1,2 and 3 stand for d, pi and P2 orbitals, respectively, to n = (2n+ 1)ttT is the Fermion Matsubara 
frequency and /x is the chemical potential. E m k a and m \ k can be calculated from e d , e p , td P and t pp . In what 
follows we express all the energies in units of t dp . 

The Dyson-Gor'kov equations for the Green's functions G\ v (k, iw n ) and the anormalous Green's functions 
Fl v (k,iu n ) are given by 

Gxu (fc, iw„) = G ( £ (fc, iw n ) + Gfl (fc, iw„)[E« (fc, iu n )Gi v (k, iw„) - £ (2) (fc, )w n )FL (fc, iw„)] , (8) 



Fl(k,iu n ) = Gi° 1 ) (-fe,-iw„)[S( 1 )(-fe,-iw„)F^(fe,ia;„) + S( 2 )(-fe,-iw„)G li ,(fe,iw„)], 



(9) 



where the self-energies due to the spin and charge fluctuations are given within the FLEX approximation as 
follows: 



£ (1) (fc,kj„) = —^2v eS (q,iil m )Gdd{k-q,kj n -iil m ), 



q.m 



T 



£ (2) (fc,kj„) = -— ^ V~sin g (q,ittm)F dd (k - q,iu n - ifi m ), 



(10) 
(11) 
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and 



with 

v cS { q ,\n m ) = u + u 2 {^Xd(q,^ m ) + \ x c M^ m ) - \{TM& m ) + t d {q^ m )}\, (12) 
v siDg (q,in m ) = u + u 2 [^Xd(q,M m ) - ^ x c d {q^ m ) - \{TM& m ) -xS(<z,i^m)}], (is) 

x-(9.A B )- 1 _ t/5a(gjinm) , ^' in ™)-i + c^c ( g,in m )' (14) 

Xd(g.i^m) = -— ^[G d(i (fc + q,iw„ +iO TO )G d( j(fe,ia;„) + F dd (fe + q, iw n + iQ, m )F dd (k, iw„)], (15) 

h.n 

Xd(q,i^ m ) = -~^rr\\[Gdd{k + q,iuj n + i£l m )G dd (k, iw„) - F dd (fe + q, iw„ + ifi m )F dd (fc, iw„)], (16) 

h.n 

where O m = 2rmrT is the Bose Matsubara frequency. 

The number of d- and p-electrons per site (respectively, rid and n p ) and the doping concentration S are 
given by 

n d = 2 lim — Ve 1 "™ T G (i d(fc,iwn), 
iVo * — ' 

n p = lim — Ve IW " T [G 2 2(A;,iw„) + G 33 (fe,iw„)], 

(5 = 5 - (n d + 2n p ), 

where the "half- filled" state corresponds to n d + 2n p = 5. 

For the purpose of calculating T c we linearize eqs.(8) and (9) with respect to F dd (k,iuj n ) or £( 2 )(fc,iw„) as 
follows: 

! 1 £ (1) (fc>«), (17) 



G dd (k,iuj n ) G d ° d (k,iu) n ) 
F^ d (k,i(j n ) = G dd (-k, -iuj n )T. {2) (-k,-iuj n )Gdd(k,iuj n ). (18) 

The normal Green's functions arc now calculated self-consistently from eqs.(10), (12), (14) and (17). The 
critical temperature for the superconductivity is determined as the temperature below which the linearized 
equation for £( 2 )(fc, iu) n ) has a non-trivial solution. The linearized equation is given by 

E (2) (fc,io)„) = -— ^2[V S in g (k - p,ioj„ - \uj m )\G dd {p,\uJm)\ 2 ]F dd -,o^ {2) {p, iw m ). (19) 

p.m 

T c is thus obtained by solving this eigenvalue problem, eq.(19); T c is the temperature where the maximum 
eigenvalue becomes unity. The superconducting order parameter has the same symmetry as Y*( 2 \k, \uj n ). 

Below T c we need to solve the nonlinear equations (from (7) to (16)) self-consistently. From the solution 
we calculate the one-particle spectra A(k,cu) and the transverse dynamical spin susceptibility \ h (q,w) by 
analytically continuing G dd (k,iu) n ) and Xd(<7>i^m) from the imaginary-axis to the real-axis. The density of 
states p(ui) and the local spin susceptibility are given as follows: 

fc 

A(k, lu) = --lm[G M (k, lu + irj) + G 22 (fe, w + ir?) + G 33 (fc, w + ir?)], (21) 

7T 

ImxH = -^^Im X - + (Q,a;), X ~ + (<7, ^) = Xd(<7, ^ + i^), (»? -> +0). (22) 
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Table 1: Parameter values for the d-p model and the calculated transition temperature T c in units of td P for 



various doping concentrations. 



e p - M 


n d 


n p 


(2 - n d )/{A - 2n p ) 


5 


T 

J- c 


-2.600 


1.314 


1.799 


1.703 


0.0890 


0.006621 


-2.575 


1.308 


1.794 


1.683 


0.1029 


0.006597 


-2.550 


1.302 


1.790 


1.663 


0.1170 


0.006475 


-2.500 


1.291 


1.782 


1.624 


0.1451 


0.006015 


-2.450 


1.280 


1.773 


1.585 


0.1738 


0.005351 


-2.400 


1.269 


1.764 


1.547 


0.2037 


0.004529 



0.0080 



0.0060 



0.0040 



0.0020 



0.0000 
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Fig. 1. Calculated values of T c for various doping concentrations. 



We define the gap functions in the following conventional way: 

uZ(k, u)) = lo — -[£W(fc,w + ir?) -E^Ob.-w -«?)]. (24) 

3 Results of Calculation 

3.1 Some details of numerical calculations. 

To study the superconducting state, we first calculate the superconducting transition temperature T c . For this 
purpose we calculate the kernel of the eigenvalue equation (19) at a fixed temperature, solving eqs. (10) and 
(12-17) for Xd(?' i^m), Xd(9: i^m) an d Gdd(k, iw„) self-consistently. Then, using the results of this calculation, 
we solve the eigenvalue problem. If eq. (19) do not have a non-trivial solution at a given temperature, we 
carry out the same calculation at a slightly decreased temperature. We repeat this routine until a error for 
T c attains 10 -6 . Below T c , we calcurate Xd(<Zji^m), Xdil^^m), G d d(k,iio n ) and Fdd(k,\uj n ) by solving the 
non-linear equations (7) to (16), self-consistently. 

All the sums involved in the self-consistent equations are carried out by using fast Fourier transforms (FFT) 
with the 64 x 64 meshes in the fc-space and the Matsubara frequency sum up to \oo n \ w 25|td p |. The solution 
is obtained by iteration until the self-consistency condition is satisfied within the following accuracy: 

| S «(*,L,„)- S «( fc ,»Ji <10 _» m=1 _ 2 . 



|E< m »(fc,iu„)| 

The analytical continuations to the real frequency axis are carried out by using Pade approximants. 

We need to take special care in evaluating p(w) and Imx(w) because of the size-effect due to the striking 
reduction of the quasi-particle damping in the superconducting state; after substituting the fc-summations 
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for Gdd(k,iuJ n ) or x^(<jr,if2 m ) with numerical integrations by using the spline interpolation method, the Pade 
approximants are calculated. 

Now, for the actual calculations, we choose the same parameter values as in I, tpp — 0.3, 6cZ £p — -1.8,1/ = 
4, in units of td p - 



3.2 Transition temperature. 

Although the present calculation has a better accuracy than the previous one, the results as shown in Table I 
and Fig. 1 arc only slightly different from those in I. The calculated value of T c and its doping concentration 
dependence are reasonable in the optimal and over-doped regimes while the observed reduction of T c with 
decreasing doping concentration in the under-doped regime is hard to reproduce. The pairing symmetry of 
the calculated state is always of d x 2_ y 2. Thus it seems appropriate to confine ourselves in what follows to the 
superconducting state of this symmetry. 
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Fig. 2. Temperature dependence of 2A k /k B T c for k = (tt,0), (3tt/4,0), (tt/2, 0) and (tt/4,0). (a) ,5 = 0.1451, (b) 
8 = 0.2037. 




Fig. 3. Anisotropy of the gap function at various temperatures, (a) 8 — 0.1451, (b) 8 = 0.2037. 



3.3 Gap function. 

We calculate the gap function as defined conventionally in the Eliashberg theory; = A(fc, A&). Fig. 2(a) 
shows the temperature dependence of 2Ak/k^T c for k = (n,0), (37r/4, 0), (tt/2,0) and (-7r/4,0) for S = 0.1451 
and Fig. 2(b) shows the corresponding results for S = 0.2037. We find that the superconducting gap develops 
more rapidly than in the BCS model. The maximum gap at k — (ir, 0) saturates to a low temperature value 
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of around 10k^T c . We next show in Figs. 3(a) and 3(b) the fe-dependence of at varying temperatures for 
S = 0.1451 and S = 0.2037, respectively. The results are roughly proportional to | cosfc^ — cosk y \. 

3.4 Density of states. 

Figs. 4(a) and 4(b) show the calculated one-particle densities of states at T — T c (dashed lines) and T — 0.0030 
(solid lines) for S = 0.1451(T C = 0.0060) and <5 = 0.2037(T C = 0.0045), respectively. It is seen in Fig. 4 
that the pseudo-gap opens at T = T c due to the antiferromagnetic spin fluctuation. The superconducting 
gaps estimated from these figures are 2A = 0.0663 for <5 = 0.1451(2A/fc B T c = 11.0) and 2A = 0.0438 for 
S = 0.2037(2A/fcBT c = 9.7). The deviation of the frequency of the minimum density of state from the Fermi 
level at T = 0.0030 for S = 0.2037 seems to indicate that a sizc-cffcct still remains in the numerical work. 
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Fig. 5. Calculated dynamical susceptibilities; upper panels show the real and imaginary parts of the irreducible 
susceptibilities and the lower panels show the final results for S = 0.1451. (a) q — (ir, n), (b) q = q inc = (ir, 0.81257r). 
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Fig. 6. Calculated dynamical susceptibilities; upper panels show the real and imaginary parts of the irreducible 
susceptibilities and the lower panels show the final results for S = 0.2037. (a) q — (it, 7t), (b) q = q inc = (it, 0.75007r). 



3.5 Dynamical susceptibility. 

We first show in Figs. 5(a) and 5(b) the w-dependence of the dynamical susceptibility in the superconducting 
state for 5 = 0.f45f for the wave vectors q = Q = (it, ir) and q = q inc , respectively, where q inc is the wave 
vector for the incommensurate peak of the susceptibility. The upper panels show the real and imaginary parts 
of the irreducible susceptibility r x s d (q : u} + irj) and the lower panels show the imaginary part of }f d {q,u) + ir;). 
The results of the same calculation for 6 = 0.2037 are shown in Fig. 6. From these figures we find that 
Imx^(q, u> + ir?) has a peak at the value of ui — 2A max corresponding roughly to the energy gap in the density 
of states shown in Fig. 2: 2A max = 0.064 for 5 = 0.1451 and 2A max = 0.049 for S = 0.2037. On the other 
hand, Imx s d (q,uj + ir/) shows a sharp resonance peak at a frequency substatially smaller than 2A max . This 
peak is considered to be specific to the superconducting state since it tends to vanish at T c as may be seen in 
Figs. 5 and 6. From the same figures we may see that the resonance peaks grow at the expense of the lower 
frequency spin fluctuations in the normal state. Also, the peak intensity seems to be stronger when the energy 
gap is larger. The resonance nature of this peak may best be seen from the following expression: 

lmXdiq > " + ^ = [l-URer d (q^ + i V W + [UlraMq,- + W ^ 

and Figs. 5 and 6. In the denominator of the right hand side of this expression the second term may be quite 
small for the energies less than the superconducting gap energy. When the first term vanishes at w = cj q , 
Inr^(<jr, u) + iry) naturally shows a sharp peak there. In Figs. 5 and 6 the resonance frequency u) q is given by 
the point where Rex^(g, u> + irj) crosses the horizontal dotted lines indicating 1/U. The peak width is given 
roughly by the magnitude of j q = Imx|(<7,Wq + iry). Figs. 7(a) and 7(b) show the calculated q-dependences 
of uj q and -f q for 5 = 0.1451 and 6 — 0.2037, respectively. It is interesting to note that ui q - q curves are 
qualitalively similar to the corresponding results for the Hubbard model. 

Next we show in Figs. 8(a) and 8(b) the q-integrated or local dynamical susceptibility Imx(^) for S = 0.1451 
and S = 0.2037, respectively. Here the solid and dashed lines show the results at T = 0.0030 and T = T c , 
respectively. 

Finally we note that the resonance peaks calculated here and the result of calculation for Im^(w) bear 



striking resemblance to the results of neutron inelastic scattering experiments on YBa 2 Cu 3 07 fl9j, Efj, El, 22 
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and YBa2Cu306,5- [P3| For more convincing comparisons, however, further experimental and theoretical studies 
are needed. 



0.06 



0.04 



3 0.02 



o.oo 



8=0.1451 
7=0.0030 



(a) 



O V 

'a 



0.06 



0.04 



3 0.02 



0.00 



8 =0.2037 
T=0.0030 



(b) 



3° * % 

° Y„ 



...22.9fies O..Q.a.Q.D°° 



(71,71/2) 



(7C.7C) 

g 



(7C/2,7C/2) 



(71,71/2) 



(7I,7C) 
<7 



(71/2,71/2) 



Fig. 7. The resonance peak positions of the dynamical susceptibility and the damping constant, (a) 8 = 0.1451, (b) 
5 = 0.2037. 
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Fig. 8. Imaginary parts of the local or q-integrated dynamical susceptibility, (a) 8 — 0.1451, (b) 8 — 0.2037. 



3.6 One electron spectral density. 

We show in Fig. 9 the calculated one-particle spectral density A(k,ui) for S = 0.1451 for five wave vectors 
around fc = (tt,0), i.e., (a) fc = (0.750tt,0), (b) fe = (0.875tt,0), (c) fe = (tt,0), (d) k = (7r,0.1257r) and (e) 
k = (ir, 0.2507r), where the solid and dashed lines show the results for T = 0.0030 and T — T c , respectively. 
Corresponding results for S — 0.2037 are shown in Fig. 10. From Figs. 9 and 10 we find a remarkable feature of 
A(k, to) in the superconducting state, consisting of a sharp coherent quasi-particle peak followed by a structure 
of dip and hump. As a matter of fact these behaviors are quite similar to those observed in the angle-resolved 
photoemission spectrum. (24|, [2^, ||^, [27) 

Next we show in Figs. 11(a) and 11(b) the dispersions of the coherent quasi-particle peaks and the humps 
together with the quasi-particle dispersion at T — T c around k = (tt,0) for 6 = 0.1451 and 6 — 0.2037, 
respectively. The dispersions are defined by the value of e(k) giving a maximum value of A(k,ui) for a 
given k. These figures show that the quasi-particle peaks are almost dispersionless while the humps become 
indistinguishable from those in the normal state as we go away from the Fermi level. This behavior is also 
consistent with the results of ARPES measurements. (27| 

In order to get better understanding of the above mentioned specific feature of the one-particle spectral 
density in the superconducting state, we show in Fig. 12 the calculated self-energy and the one-particle spectral 
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Fig. 9. One-electron spectral density around k — (ir,0) for 8 — 0.1451. 
Fig. 10. One-electron spectral density around k = (tv, 0) for 8 — 0.2037. 
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Fig. 11 Energy dispersions of the coherent quasi-particle peak and hump in the superconducting state (T = 0.0030). 
(a) 6 = 0.1451, (b) S = 0.2037. 



density of the d-state, since the spectral density is dominated by this contribution. We have 

1 1 



G dd (k, uj + in) G { °] (k, oj + it]) 



£(fc,w + irj), (27) 



•dd ' 

A d (k,u)) = --lmG dd {k,u) + ir]). (28) 



E(fc,w + ifj) = Z w {k,LJ + iri)-G^{k,-uj + iT])[rt 2 \k,u + i7])} 2 , (29) 

EWflfe.w + ty). (30) 



G N (fc,^ + i?7) G^(fe,a; + i7y) 
A quasi-particle peak is expected at the frequency which satisfies 

1 



RcZ(k,uj + ir)) = 0, (31) 



provided ImE(fe,w + iry) is small at the same time. As is seen in Fig. 12 the superconducting energy gap 
gives almost vanishing values of ImS(fe,w + irf) in a fairly wide frequency range around k = (-7T, 0) and sharp 
quasi-particle peaks are expected in this region of fc-space. As the frequency increases beyond this region 
|Im£(fc, u> + ir))\ increases first rapidly and then gradually. At the same time the value of [G d °J (k, lo + ir/)] -1 — 
RcS(fc,o; + irf) increases rapidly first to form a dip structure in A d (k,u>) and then the ^-dependence of this 
value reflect itself on the spectrum of A d (k,ui) as is expected from eqs. (27, 28) and is actually seen in Fig. 
12. 



4 Conclusion and Discussion 

On the basis of the d-p model we have studied the unconventional superconducting state induced by the 
antiferromagnetic spin fluctuations in the optimal and over-doped regimes of high-T c cuprates. Calculations 
were performed for the superconducting oder parameter of d x 2_ y 2 symmetry, the dynamical susceptibility 
and the one-electron spectral density, by using the FLEX approximation, i.e., the strong coupling theory 
with the self-consistently renormalized random phase approximation (RRPA) for spin fluctuations. The order 
parameter or the gap function was found to develop more rapidly than in the BCS model, and the ratio of 
the low temperature maximum gap energy to k^T c was about 10. The calculated spin fluctuation spectrum 
below T c showed a resonance peak around q = (tt, it) whose intensity was especially strong around the optimal 
doping concentration and its behavior bears striking resemblance to the peak at 41 meV observed in the 
neutron scattering experiment on YBa 2 Cu 3 07. The calculated one-particle spectral density showed a sharp 



10 




Fig. 12. Self-energy of the d-component of the Green's function and its one-particle spectral density at three points in 
the fe-space around k = (tt, 0). 



quasi-particle peak around k = (w, 0) followed by a structure with a dip and a hump. These results were quite 
similar to those observed in the ARPES spectra. 

In view of these results we may conclude that the present calculation gives an additional support to the 
spin fluctuation mechanism for the high-T c cuprates, especially in the optimal and over-doped concentration 
regimes. 

Desired improvements of the present approach include the consideration of the vertex corrections neglected 
in the FLEX approximation. Although the effect was suggested to be insignificant or rather favorable for 
the superconductivity, |2a] a self-consistent calculation with the vertex corrections is desired. The effect of 
impurity potential is expected to be more significant in the superconducting state than in the normal state. 
For example, the observed linear temperature dependences of 1/T\ and the specific heat far below T c were 
successfully interpreted in terms of the additional density of states around the Fermi level caused by the 
impurity effect. [p9[ This effect will also have significant influence on the resonance peak in the dynamical 
susceptibility and on the above-discussed structures in the one-electron spectral density. 

Of course there are many problems still to be clarified, particularly in the under-doped regime. The most 
controversial among them may be the pseudo-gap phenomena in the under-doped regime observed in NMR, |30| 



ARPES, plL |32|, 33, B4| |35| and other measurements. In the present approach we treat the spin fluctuations 
in the spirit of the self-consistent renormalization (SCR) theory which applies around the antiferromagnetic 
instability. However, the FLEX approximation without the vertex corrections may not be sufficient particularly 
in dealing with the low frequency phenomena. Also, it seems possible that the ground state in the under- 
doped regime has an antiferromagnetic long range order extending to the region rather far from the instability 
point. If such circumstances are really the case it seems necessary to develop an approximation beyond the 
conventional SCR, taking account of substantially larger effects of non-linear mode-mode couplings. 
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